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As a most promising candidate for quantum theory of the gravity, the superstring 
theory has attracted many researchers including cosmologists. It is expected that 
the cosmological initial singularity is avoided within the context of the superstring 
theory. Indeed, Antoniadis et. al. found an interesting example of the the non- 
singular cosmological solutions by considering the string one-loop correction. Here, 
we will discuss the stability of this model and find a new kind of instability in 
the graviton modes. We also argue this instability will persist even in the non- 
perturbative regime. The instability we have found might provide a mechanism to 
produce the primordial black holes. 



1 Introduction 

As is well known, the singularity is unavoidable consequence of the general rela- 
tivity. The notorious cosmological singularity is nothing but the manifestation 
of this general theorem. Usually, we expect that the quantum gravity works 
at the Planck scale and the singularity disappear. To date, the superstring 
theory is the most promising candidate for the consistent theory of quantum 
gravity. Hence, it is interesting to seek the non-singular cosmological solutions 
in the low energy effective action of the superstring theory. From this point of 
view, Antoniadis et. al.ld studied the one- loop corrected heterotic superstring 
effective action and found a particular interesting class of solutions which avoid 
the initial singularity. 

The purpose of this work is to show that these solutions are unstable due 
to the high frequency gravitational wave. It means that the background geom- 
etry is not viable in the most cases. When we allow the fine tuning, however, 
the instability leads to an interesting consequence. If, at the arbitrary initial 
time, the amplitude is sufficiently small, then the small scale high amplitude 
gravitational waves will be produced by the instability without destroying the 
background geometry. In the deceleration phase, they will collapse and become 
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the primordial black holes. To support the above picture, we will discuss the 
non-perturbative effects phenomenologically and demonstrate that the insta- 
bility will persist even in the non-perturbative regime. 

2 Non-singular Cosmological Solution 

From the non-linear sigma model approach, we can calculate the (3 function 
of the coupling function. The Weyl invariance requires the vanishing of the (3 
function, which gives the equations of motion for the massless modes. From 
the equations of motion, the effective action can be easily read off. When we 
fix the world surface as the sphere, we obtain the tree level effective action. In 
actual calculation, we usually expand the effective action with respect to the 
string tension which controls the field theoretical loop expansion. The string 
loop correction is different from this. The string loop expansion corresponds 
to the path-integral quantization on the higher genus manifold. In particular, 
the one-loop string correction can be obtained from the calculation of the f3 
function on the torus. It is this kind of correction that we want to investigate 
in detail. 

The jaction for the heterotic superstring with orbifold compactification is 
given by El 

S = J V=9~{\r ~ \{D<S>? - \{Dof + ^[A ie * _ X 2 i{a)]Rl B ) (1) 

where R, $ and a are the scalar curvature, dilaton, and the modulus field, 
respectively. The Gauss-Bonnet curvature is given by 

Rgb — R^^Rpvpo — ^R^Rfiv + R 2 (2) 

and the coupling function can be written by using the Dedekind rj function: 

$(a) = -log[2e° V (ie°)] (3) 

For the existence of the non-singular cosmological solutions, the modulus field 
a is relevant. Therefore, from now on, we will ignore the dilaton part aLthc 
action. Thus, we shall examine the modulus part of the effective actionem 

S = J d^x^—g^-R - \{D<PY - ^H(4>)R% B } . (4) 

Under the homogeneous and isotropic ansatz 

ds 2 = -dt 2 + a 2 {t)5 ij dx i dx j , (5) 
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we obtain the equations of motion: 

-zh* = - 1 -4?-^hH (6) 

-[2H + 3H 2 ] = ^-^[H 2 i+2Hi(H + H 2 )} , (7) 

where H = a/a = da/dt/a. 

It is convenient to define the effective energy density and pressure as p e s = 
-G° and p cS = G°/3. Then, 

r Peg + Pes 2H 

Here, we would like to stress that this effective adiabatic index is that of the 
background effective matter. It is not related to the adiabatic index of the 
pcrturbative matter. 

Utilizing the asymptotic form 

O~sign(0)|el*l, (9) 

it is possible to obtain the asymptotic solutions by imposing the ansatz 

H = Lo 1 \tf (10) 
= 0o+w 2 log|i| . (11) 

In the asymptotic future, we expect the Friedman-Robertson- Walker universe, 
hence assume that the Einstein part balances with the ordinary energy mo- 
mentum in the future. The asymptotic solution in t — > oo becomes 

a ~ , (12) 

"'Vl- < 13 > 

This is the Friedman-Robertson- Walker solution with stiff matter. For the 
later purpose, we need to calculate the effective adiabatic index. From the 
asymptotic solutions, it is easy to calculate T as 

r = 2,<^oo. (14) 

In the asymptotic past, t — > — oo, in order to obtain the non-singular 
solution, the Gauss-Bonnet part must balance with the ordinary energy mo- 
mentum. Then, we obtain the super-inflationary phase 
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This type of solutions appears in the pre-big-bang scenario. The serious prob- 
lem there is the graceful exit problem. Interestingly, one-loop corrected action 
gives a solution to this problem. The effective adiabatic index becomes 




t 



(16) 



Numerical calculation tells us that the above asymptotic solutions are smoothly 
connected. This is an example of the graceful exit in the pre-big-bang sce- 
nario. tl 

Apparently, in the Gauss-Bonnet dominated phase, the energy condition 
is violated strongly. This fact itself is expected generally to obtain the non- 
singular cosmological solutions. 

3 Stability analysis 

Intuitively, it is expected to have the instability in the background which breaks 
the energy condition. In case of the perfect fluid, the negative pressure causes 
the instability and the inhomogeneity is enhanced catastrophically. However, 
what kind of instability can we imagine in this field theoretical context? Does it 
happen in the scalar perturbation, vector perturbation, or tensor perturbation? 
Unfortunately, we do not have any answer without calculation. Indeed, we 
performed numerical calculations and found no essential instability in the scalar 
and vector perturbation. It is the tensor perturbation mode that shows the 
instability. Here is something interesting. The isotropic and homogeneous 
background drives the instability of the gravitational wave modes. This counter 
intuitive result prevents us to use the fluid analogy. After all, the analogy is 
nothing but analogy. Another interesting feature of this instability is its scale 
dependence. On the contrary to the Jeans instability, the instability becomes 
strong in the high frequency graviton modes. Hence, the time scale of the 
instability could be arbitrarily small as far as the cut off scale is not introduced. 

To illustrate the instability explicitly, it is convenient to see the Hamilto- 
nian of the system. Consider the tensor perturbation 



ds 2 = -dt 2 + a 2 ((5.y + h ij )dx i dx j (17) 

where h 1 ^ = h\ = 0. By using the formula listed in the Appendix, we obtain 
the action for the tensor perturbation as 
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-£[-a 3 Hh,jh lj + Aa 3 (H + H 2 )Hh i3 h 1: >}} (18) 
= \J d 4 xa 3 [(l - ^Hbhiih* - (1 - ^)^V/^W] (19) 

The last expression is obatined by using the background eqations. From the 
background equation, we also get 

l-^ = (l-^)(-f -5) 

= (l-^ff«)(3r-5) (20) 

Notice that a = 1 — > 0. Finally, we have obtained the action for the 

graviton modes: 

S = IJ tPxcPaihijh*' - (3r - 5)^V/iyV/i y ] (21) 

It is straightforward to deduce the Hamiltonian of the system in the following 
form: 

H = [ d 4 x[^^- + (3r - 5)aaVh ij Vh ij ] , (22) 
J a 6 a 

where 7r y = a 3 ah^ /4. 

If the energy condition is violated, i.e. r < or T < |, the system is 
unstable. Instability becomes strong as the wavelength becomes short. 

In case of one-loop model, we have T = 2 in the Friedmann-Robertson- 
Walker phase, so the background geometry is stable as is expected. On the 
other hand, it turns out that the Gauss-Bonnet phase is unstable due to the 
strong violation of the energy condition, T ~ — \t\. Strictly speaking, this insta- 
bility immediately does not imply the breakdown of the background geometry. 
The point is that the fine tuning of the graviton amplitude is required to avoid 
the breakdown of the background geometry. It is possible Xo give the criterion 
for the breakdown of the background geometry explicitly. □ 

4 Discussion 

First of all, it should be noted that the appearance of the Gauss-Bonnet term 
is generic features of the superstring cosmology. In case of the one-loop correc- 
tion, the function form of the coupling between the modulus and the Gauss- 
Bonnet term is universal. Hence, our result can be regarded as the general 
feature at this level. However, our result suggest that the non-perturbative 
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effects should be taken into account if superstring theory is viable for solving 
the singularity problem. Assuming the Damour-Polyakov type universality, we 
can obtain the non-perturbative effective action 

S = J d 4 x^g{R - (d$) 2 + B{$)R 2 GB + ■■■} (23) 

where 

B(*) - e* + c + cie"* + c 2 e^ + ■■■ . (24) 

Here, the function B($) contains the correction from the string higher loops. 
The coefficients Cj should take the values so that the initial singularity will be 
avoided. The energy condition is necessarily violated in this case again. Now, 
we shall show, even in this case, the instability we have found will appear. 
More generally, we can take into account other massless fields in the model: 



<£x^—g\R £(V0 4 ) 2 - \c{fa)R% B \ (25) 



where G(fa) is now a function of massless fields fa. The action for the tensor 
perturbation again takes the following form 

Sperturb = \ j #xa\\ - ^ HG{fa)) [hij ti j - (3T - 5)^WhijWh ij ] . (26) 

This is essentially the same as the previous form. So, we have proved the 
existence of the instability. 

Next, we will argue that the fine tuning problem is not avoidable in the non- 
singular cosmology in this model. Let us assume the non-singular cosmological 
solution. Then, we have the eternal past and the eternal future. Of course, 
there exists a period where the energy condition is violated in order to have 
the non-singular cosmological solution. However, it is not allowed to have the 
infinite period of violating the energy condition, otherwise we will encounter 
the fine tuning problem of the graviton amplitudes. This implies that the 
effective adiabatic index T is positive , more precisely T > 5/3 in most of the 
time. From the definition T = —2H/3H 2 , we obtain 

S<ffk» - 5 r « < 27 » 

As r(t) is almost always positive, 1/H(t) is almost monotonic function, hence 
1/H(t) becomes zero at a certain time to. This means 

lim H(t) = oo (28) 

t— 
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namely wc have encountered the singularity. Thus, as far as the non-singularity 
and the avoidance of the fine tuning is assumed, we must have the singularity 
somewhere. This is the contradiction. Hence, if we want to have the non- 
singular cosmological solutions in this context, the fine tuning problem can not 
be avoidable. Hence, this kind of non-perturbative correction can not change 
the situation. It should be noted that the spatial curvature might change the 
arguments above. □ 

An interesting consequence of our result is that the instability is stronger 
in the small scale. Hence, the instability found by us can provide a mechanism 
to generate the primordial black holes. We will investigate this possibility in 
the future. 
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Appendix 

Here, we will list the necessary formula to calculate the perturbative action: 

(H + H 2 )[~h ik h ik - \ti k h lk - Hh lk h lk ] 

2H 2 [——^h lk ' m hik,m — —,'h lk 'hik — 2Hh lk hi k ] 
4:a z 4 

^h ik h ik + Hh ik hik + H 2 h ik h ik 

_J_iik,m L 
n o ,tj IL ik,rn 

2a z 

\v 2 h lk V 2 h lk - ^h lk V 2 h lk + H 2 h lk h lk 
a 4 ar 

3(H + H 2 )[-^h ik h lk - \h lk h lk - Hh tk h lk ] 
(H + iH 2 )[-^h lk h lk - \h lk h lk - 3Hh ik h lk ~ ^h jk>m h jk ' m ] 
h tk h lk + ^V 2 h ok V 2 h^ k + ^H 2 h lk h lk 

-^2^kV 2 h jk + ^Hh lk h tk - ^h jk V 2 hi k ] 
6(H + 2H 2 )[-~h ik h ik - h lk h lk - ±Hh lk h lk - ^hft^] . 



„(o)cy p (2)m 

K ij K kl 

R (i)0i „(i)fci 

K kl K 0j 

R (i)y p(i)fc* 

d(0)0 p(2)0 
T>(0)i r>(2)j 

Kj Hi 
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Here, we have omitted the spatially total derivative terms, because they do 
not contribute the final answer. 



References 

1. I. Antoniadis, J. Rizos and K. Tamvakis, Nucl. Phys. B 415, 497 (1994). 

2. J. Rizos and K. Tamvakis, Phys. Lett. B 326, 57 (1994). 

3. M. Gasperini and G. Veneziano, Astropart. Phys. 1 (1993) 317. 

4. S. Kawai, M. Sakagami and J. Soda, Instability of 1-loop superstring 
cosmology, |gr-qc/9802033 . 

5. R. Easther and K. Maeda, Phys. Rev. D 54, 7252 (1996). 



8 



